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Using a simple lumped-circuit model, we numerically study the dependence of 
the voltage gain and noise on the amplifier's parameters. Linear, quasi-linear, and 
nonlinear regimes are studied. We have shown that the voltage gain of the amplifier 
cannot exceed a characteristic critical value, which decreases with the increase of the 
input power. We have also shown that the spectrum of the voltage gain depends 
significantly on the level of the Johnson noise generated by the SQUID resistors. 



I. INTRODUCTION 



A microstrip-SQUID (superconducting quantum interference device) amplifier (MSA) has 
been designed as a low noise radiofrequency amplifier, which is able to operate above 100 
MHz [l|. The MSA has been studied theoretically in many publications [l 9 ] . However, 
a consistent theoretical model of MSA has not yet been developed. The circuit diagram 
of lumped model of MSA is presented in Fig. [TJ This MSA consists of a linear input 
circuit coupled to the direct current (dc) SQUID via the mutual inductance, M. Note 
that the isolated dc SQUID is a nonlinear circuit, while the isolated microstrip is a linear 
circuit. The total system consisting of the SQUID and the input circuit is a nonlinear 
one. Consequently, it remains no trivial task to predict the performance of the MSA during 
the design process. Therefore analytical investigation and extensive numerical modeling, 
simulation, and optimization of the MSA are required before creating the device. Generally, 
the solutions for the MSA model must be obtained by solving analytically or numerically 



the system of nonlinear ordinary differential equations written for the SQUID coup 



the input circuit. The dynamics of a bare SQUID was investigated numerically in 1CH12] . 



ed to 



The objective of our paper is to study numerically the dynamics of exact nonlinear equa- 
tions describing the MSA, and to compute the voltage gain, G(f) = \V(f)/Vi\, of the MSA, 
where V(f) is a Fourier harmonic of the output voltage on the SQUID, and V% is the ampli- 
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Input circuit 

FIG. 1: The equivalent scheme of the microstrip-SQUID amplifier. Vi is the amplitude of the 
input voltage and R\ is the resistance of the voltage source. L\ describes the stray inductance, 
and (if necessary) the pick-up coil inductance, G\ is the coupling capacitance, L is the inductance 
of the input coil for the SQUID, / and if are, respectively, the bias current and the flux for the 
SQUID, M is the mutual inductance between the input coil, L, and the SQUID loop. Lj is the 
inductance of the SQUID loop, I\, Ic, and Ij are the currents in the input part of the circuit, and 
J is the current circulating around the SQUID loop, Cj and Rj are, respectively, the capacitance 
and resistance of each Josephson junction. 

tude of the input voltage on the microstrip at the same frequency, /. We analyze numerically 
both linear and nonlinear regimes of amplification. A linear regime means that the following 
linear dependence exists: |V(/)| = G(f)\Vi\, where the gain G(f) is independent of the Vi. 
We also simulate the output spectral density of voltage Johnson noise, originated in shunt- 
ing resistors of the SQUID and the resistor, R±, in the input circuit, and calculate the noise 
temperature. 



II. INPUT CIRCUIT 

Consider the isolated linear input circuit (M = 0). The forward impedance of the input 
circuit is 

Zi = ^ = (Wi/Li + — ^— + Ri) (l " (2tt/) 2 LC) + 2mfL. (1) 

In Fig. [2] we plot the amplitude of the current, Jj, in the input coil with inductance L, 
when M — 0. One can see from Fig. [2] that as C\ decreases, the maximum shifts to higher 



3 



0.20 



0.15 



% 

^ 0.10 

o 

0.05 



0.00 







C=0.5 pF 




C=2.2 pF 









100 200 300 400 500 600 700 800 900 1000 

Frequency (MHz) 

FIG. 2: Current, Ii, in the input coil for an isolated input circuit (M = 0), and for three values 
of the coupling capacitance, C\. R\ = 50 Ohm, V% = 1 /iV, L\ = 1 nH, C = 4.4 pF. /o = 
l/(27r\/ LC) = 700 MHz is the resonant frequency of the input resonator. 



frequencies, and the width of the peak decreases. One can use this latter property to create 
a narrow-bandwidth amplifier. The frequency corresponding to the maximum is always less 
than the resonant frequency, / , of the input resonator. 



III. EQUATIONS OF MOTION 

The differential equations of motion for the SQUID are: Q 

(foCjSx + h = I- - J - I sin 5i + I nl , 
Kj I 

(PoCj5 2 + ^-5 2 = - + J - I sin 5 2 + I n2 , (2) 
Kj z 

<Po(S 1 -6 2 ) = ip + LjJ + MIi. 

Here the dot above 8\ and 5 2 indicates time differentiation; ip Q = %/{2e) is the reduced flux 
quantum; h is Planck's constant; e is the electron charge; 5\ and S 2 are the phase differences 
in the Josephson junctions in the SQUID; Iq is the Josephson junction critical current; l n \ 
and I n2 describe the noise current (Johnson noise) originating in the shunt resistors, Rj. The 
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output voltage, V, can be expressed in terms of the Josephson junction phase differences 

v = f (4 + 4). 

In order to determine Jj in the third equation in (T2J), we have to add the differential equations 
for the input circuit. The total system of eight first-order differential equations can be written 
in the following form: 



Here 



Pi 

P2 

Qi 



Pi, 



P2, 
1 



R 



-P2 



2 sin - — cos h I nl + I n2 

2(p 2(fo 



(x 2 - ip — Mli) - 2 sin cos + I nl - I n2 



2<p 2<Pc 



V, cos(27r/t) + V n -%- hRx - % 



Qc 

l = 



h - h, 

1 (Qc 



M 

Ln V C ' T/ 2 



^1 = ^0(^1 + ^2); ^2 = ^0(^1-^2); 



a 



1 - 



M 2 
LL, 



(3) 



Qi is the charge on the capacitor C\\ Qc is the charge on the capacitor C; / is the frequency 
of the external voltage; and V n = I n Ri is the noise voltage on the the resistor Ri, I n is the 
noise current through R\. 

The input circuit is coupled to the SQUID through the term, MIi/(CjLj), in the fourth 
equation in ([3]), and the SQUID is coupled to the input circuit by the effective coupling 
constant, 7 = M/(LLja), in the last equation in ([3]). Since the effective coupling constant, 
7, is proportional to 1/a, the effective coupling can be increased by decreasing a. 



IV. VOLTAGE GAIN 



First we choose the optimal working point of our amplifier which is defined by the value 
of <p, provided the other parameters are given. In Fig. [3] the time average, V , of the output 
voltage and its derivative (transfer function) are plotted as a function of ip. The maximum 
of the transfer function occurs in the vicinity of (p/(fo — 0-3 which we choose as our working 
point. 
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FIG. 3: Average over time of the output voltage, Vo, (left) and its derivative (transfer function) 
(right). The working point for 92 /Vo = 0.3 is marked by a filled blue circle. Cj = 0.2 pF, Lj = 0.45 
nH, Rj = 20 Ohm, I = 1.99I , I = 2 fj,A, a = 0.001, C x = 0.5 pF, V- = 0, I nl = I n2 = 0. The 
other parameters are the same as in Fig. [2j 
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FIG. 4: (a) Voltage gain of the amplifier for three values of the coupling capacitance, Ci, and for 
a = 0.001. (b) Voltage gain for three values of a, and for C\ = 0.5 pF. Vi = 0.001 fiV in both (a) 
and (b). The other parameters are the same as in Figs. [2] and El 



In Fig. @|a) we plot the voltage gain for three values of the coupling capacitance, G\. 
By comparison of Fig. H)^a) with Fig. |2j one can conclude that the amplifier gain is mostly 
defined by the parameters of the input circuit. The gain of 55 dB for C\ = 0.5 pF is 
the combined result of amplification by the input circuit and by the SQUID due to the 
strong interaction between them when a is small and positive. Note that the dimensionless 
parameter, a, contains both the parameters of the SQUID (Lj) and the input circuit (L), 
as well as the coupling inductance, M. In Fig. H^b) we plot the gain for three values of 
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FIG. 5: Voltage gain for three values of input voltage amplitude Vi. a = 0.001, C\ = 0.5 pF, the 
other parameters are the same as in Figs. [2] and [3l 



a: a = 0.001, 0.01, and 0.1. The gain decreases from 55 dB to 27 dB as a increases from 
0.001 to 0.01. The value of coupling inductance, M, changes, respectively, from 2.2982 nH 
to 2.2878 nH, that is, by only 0.44 percent. Therefore, the possibility of obtaining the large 
gain is limited to a very small region of the parameters. For this purpose it is desirable 
to have a tunable coupling inductance, M, or a tunable input circuit inductance, L, or a 
tunable SQUID inductance, Lj. 

Consider the situation in which a is negative, that is M 2 > LLj. In order to understand 
the dynamics in this regime, we differentiate the last equation in ([3]) and use the 4th and 
7th equations for Qc and p2- We obtain an equation for the oscillations of the input current, 
Ii, with an external force and with the eigenfrequency, uq, where 

„ 1/1 M 2 N 

a 



Wn = - 



LC + L 2 LC. } 



Negative a corresponds to negative reai 
can drive the resonator into instability 



part of the input impedance 0, [jjj . This regime 
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V. NONLINEARITY 

When the input voltage amplitude or gain becomes sufficiently large, the nonlinear effects 
in the SQUID become important. In the nonlinear regime, the nonlinear effects decrease the 
output voltage of the SQUID, thus decreasing the gain of the amplifier in comparison with 
the linear regime. It is reasonable to assume that the MSA is in the nonlinear regime when 
the output voltage becomes comparable with the SQUID's own average output voltage, V . 
It is convenient to define the maximum gain 



G, 



V 
V 



(4) 



When the amplitude, |V(/)|, of the amplified output Fourier harmonic approaches |Vo|, the 
gain should decrease due to nonlinear effects in the SQUID. If, for example, the amplitude of 
the input signal is Vi = 1 /iV and V = 7 /iV (see left side in Fig.[3]for <p/(po — 0.3), according 
to Eq. (p|J the gain is limited by the value G max = 7. In order to obtain a gain of G = 55 dB 
in Figs. Hta) and (b), we set the amplitude of the input signal to be Vi = 0.001 fiV. For this 
input signal expressed in decibels [G max (dB) — » 20 log 10 G max ] , we have G max = 76.9 dB, so 
that G(f) < Gmax- For Vi = 0.01 /xV, we have G max = 56.9 dB, and for Vi = 0.1 /iV one 
obtains G max = 36.9 dB. In Fig. Owe plot the gain, G(f), for three values of V*. One can 
see from the figure that G(f) is less than G max for all frequencies, /. (The values of G max 
for each V, are indicated in the figures by the horizontal lines.) In Fig. 0(a) we plot G as 
a function of V.. As follows from the figure, the condition G < G max is also satisfied for 
all input voltage amplitudes, Vi. As the amplitude, Vi, of the input signal increases, both 
G m ax in Eq. (j3J) and the gain, G, decrease. This nonlinear effect cannot be obtained from a 



linearized theory, like that based on an effective impedance of the amplifier 16] . In Fig. [6(b) 
the amplifier is in the linear regime for C\ = 2.2 pF and C\ = 10 pF because the gain is 
sufficiently small, G <C G max . 



VI. NOISE 

We assume that the Johnson noise voltages across the resistors dominate all other sources 
of noise in the amplifier. White Gaussian noise in each of three resistors was modeled by a 
train of rectangular pulses following each other without interruption. The current amplitude, 
I n , of each pulse is random with zero average and the following variance: 

II 2 ) = 2 M. 
Kn) RAt 
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FIG. 6: (a) Voltage gain, G, and maximum gain, G max ^ a function of the input voltage amplitude, 
Vi. In maximum G = 58.7 dB (black curve), (b) Voltage gain, G, in linear regime as a function of 
the input voltage amplitude, Vi. a = 0.001 and the other parameters are the same as in Figs. [2] 
and [3l 

Here the brackets indicate an average over different realizations; kg is the Botzmann con- 
stant; T is the temperature of the MSA; At is the duration of each pulse, which is constant 
in our simulations; R is the corresponding resistance: R = Rj for the shunting resistors in 
the SQUID and R = Ri for the resistor R 1 in the input circuit. In Fig. [7] we plot the output 
voltage spectral density, Sy(f)/S v (where S v = (poIoRj), when no input voltage is applied 
to the input circuit, Vi = 0. Sy is defined as: 

S v = lim hv{f)\ 2 ), V(f) = f V^e^dt, 

where ti is the time of integration of the output signal and the angular brackets, (), indicate 
an average over different realizations. 

One can use the voltage spectral density, Sv(f), to calculate the noise temperature of 
the amplifier, 7\, using the following equation: 

4MT 1 +T) J R 1 |G(/)| 2 = M/). 

In Fig. Owe plot T\/T q , as a function of frequency, /, for two different scales, where T q = 
hf /ks is the quantum temperature. We use the noise spectral density, Sv(f), from Fig.[7]and 
the gain from Fig. H](a). At the minimum, for C\ = 0.5 pF (red line) the noise temperature 
is negative. 

We now show that this negative noise temperature is related to the method of its calcula- 
tion. The reason for the negative noise temperatures is in the different methods of calculation 
of the noise spectral density and the gain of the amplifier. There are three resistors that act 
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FIG. 7: Output voltage spectral density Sy/Sy, where S v = ip$IoRj, as a function of frequency, /, 
for three values of C\, a = 0.001. For comparison, we also include the plot of Sv/S v for bare SQUID 
(when M = 0) with white noise spectral density (because the frequency, /, is much smaller than the 
Josephson frequency). The time of integration is U = 10 4 io> where to = <A)/(^o-R/) = 8.23 x 10~ 12 
s; the number of realizations is N = 100; T = 100 mK; At = O.lto; Vi = 0; the other parameters 
are the same as in Figs. [2] and El 



as sources of noise: two resistors, Rj, in the SQUID and the resistor, R\, in the input circuit. 
All three resistances contribute to the noise spectral density, SV. The noise current in the 
SQUID, besides contributing to noise spectral density, changes the SQUID's parameters, 
including the frequency at which the gain is maximum. At one moment the SQUID is tuned 
to one frequency and at the next moment it is tuned to another frequency. The noise voltage 
generated in the resistor, Ri, is amplified by the detuned SQUID. On the other hand, the 
signal is amplified by the SQUID tuned to a definite frequency because there are no noise 
currents in the resistors, Rj, of the SQUID. Consequently the amplification of the signal is 
greater than the amplification of the noise. Besides, the noisy SQUID shifts the frequency 
at which the amplification is maximum. In Fig. H^a) the maximum is at f m ax = 450 MHz, 
while in Fig. [7] the maximum is at 510 MHz. 

By this argument, we calculated the gain (Fig. [H]) and the noise temperature (Fig. [TO]) 
with noise on the resistors, Rj, and with the input signal, Vi ^ 0. The calculated gain 
in Fig. [H]is qualitatively similar to that measured experimentally in Ref. 17|. If the time 
of integration is sufficiently long, the contribution of the noise to the gain is minimized, 
so only the contribution from the amplified signal remains (see Fig. [TT]) . In this situation 
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FIG. 8: Noise temperature T\/T q for three values of G\ in two different scales. The other parameters 
are the same as in Fig. E(a) and Fig. [7J 



20 



CD 

T3 



10 



'ra 5 
CD 

03 

cn n 

CD 

-t— ' 

o 

> 

-5 
-10 



1 1 1 I I I I I I i i i 



C=10pF 




C,=2.2 pF / 






fjc=0.5 pF V 



200 300 400 500 600 700 

Frequency (MHz) 

FIG. 9: Gain of the MSA with the noisy SQUID. U = 10%; V { = 0.1 /iV; N 
T = 100 mK. The other parameters are the same as in Fig. Eta). 



10; a = 0.001; 



both the voltage spectral density and the gain are calculated for the same system with the 
noisy (detuned) SQUID. Note that a significant asymmetry of the spectrum for C\ = 2.2 pF 
and 10 pF in Figs. 8 and 10 appears because we plot the ratio Ti/T qi where the quantum 
temperature T q is proportional to the frequency /. 

We demonstrated two methods of calculating gain and noise temperature of the MSA. The 
noise on the shunting resistors of the SQUID reduces gain of the amplifier if one compares 
Figs. ID^a) (no noise) with Fig. [9] (with noise). The reduction is large for small capacitance, 
G\ = 0.5 pF, while the gain for G\ = 10 pF is mostly not affected by the noise in the 
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FIG. 10: Noise temperature of the MSA with the noisy SQUID in two different scales. The 
parameters are the same as in Fig. [9j 




FIG. 11: Gain of the MSA with the noisy SQUID calculated using different integration time, t%. 
The contribution to the gain from the noise on the SQUID decreases as the integration time, ij, 
increases. Forti/to > 10 6 the gain is mostly independent of which indicates that only contribution 
from the signal remains. C\ = 0.5 pF; the other parameters are the same as in Fig. [9l 



SQUID. The gain calculated in the previous sections of this paper is actually the gain at 
zero temperature. 

In summary, we have simulated the dynamics of the microstrip-SQUID amplifier in both 
the linear and nonlinear regimes and studied the dependence of the voltage gain and noise 
on the parameters of the amplifier. We have shown that the voltage gain cannot exceed the 
critical value G max given by the formula fll]). This value is inversely proportional to the input 
voltage. It is shown that the gain decreases as the device temperature increases. Finally, we 
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have shown that the spectrum of the voltage gain depends significantly on the level of the 
Johnson noise in the SQUID resistors. This effect must be taken into account for correct 
calculation of the amplifier noise temperature. The next important step should be the 
optimization of the gain and noise temperature with respect to the amplifier's parameters. 
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